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When a symmetry gets spontaneously broken in a phase transition, topological defects 
are typically formed. The theoretical picture of how this happens in a breakdown of a 
global symmetry, the Kibble-Zurek mechanism, is well established and has been tested in 
various condensed matter experiments. However, from the viewpoint of particle physics 
and cosmology, gauge field theories are more relevant than global theories. In recent 
years, there have been significant advances in the theory of defect formation in gauge 
field theories, which make precise predictions possible, and in experimental techniques 
that can be used to test these predictions in superconductor experiments. This opens 
up the possibility of carrying out relatively simple and controlled experiments, in which 
the non-equilibrium phase transition dynamics of gauge field theories can be studied. 
This will have a significant impact on our understanding of phase transitions in the early 
universe and in heavy ion collider experiments. In this paper, I review the current status 
of the theory and the experiments in which it can be tested. 

1. Introduction 

The formation of topological defects in phase transitions is a very generic phe- 
nomenon in physics. It can be studied experimentally in different condensed matter 
systems, but it is also believed to have happened during the early evolution of the 
universe, ^'^ and the produced topological defects may have observable consequences 
to the properties of the universe today. Defect formation also has many similarities 
with the non-equilibrium phenomena that take place in heavy ion collisions. "^'^ 

In this article, however, we discuss defect formation in field theories as a phys- 
ical phenomenon in its own right, as a unique way of obtaining information about 
the non-equilibrium dynamics of field theories. Theoretical calculations of phase 
transition dynamics are only possible in practice if one is willing to make many 
simplifying approximations and assumptions, and it is not at all clear which of 
these are actually justified. In this respect, it is useful that quantum field theo- 
ries are applicable very generally in physics. Very similar field theories describe a 
wide range of systems from superconductors to interactions of elementary particles. 
In particular, the problems faced in theoretical calculations are identical in these 
systems. This raises the hope of testing some theoretical predictions in condensed 
matter experiments,^ in order to find out which approximations are reliable and 
how one should treat the field theories of particle physics theoretically. 
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To a certain extent, this hope has aheady been reaUzed in equihbrium field the- 
ory. Critical phenomena can be studied experimentally in many condensed matter 
systems, and they form the basis for the non-perturbative picture of quantum field 
theories.® In particular, the concept of renormalization group, which is based on the 
critical behaviour of field theories near second order transitions points, is crucial for 
our understanding of particle physics. 

In the case of non-equilibrium dynamics, both experiments and theoretical cal- 
culations are significantly more difficult, and therefore less progress has been made. 
In particular, it is very difficult to study directly the non-equilibrium dynamics dur- 
ing a phase transition in any experiment. Topological defects play a special role in 
this respect, because they survive for a long time after the phase transition, when 
the system has otherwise equilibrated. Their number density and spatial distribu- 
tion carry information about the dynamics of the phase transition, and they can be 
compared to predictions of theoretical calculations. 

In many ways, this approach is very similar to particle cosmology, where we hope 
to extract some information about particle pliysic;s from the relics of the Big Bang, 
and to particle accelerator experiments, in which we cannot observe the collision 
itself, but by studying its products we can reconstruct it. The difference is that 
superconductors and other condensed matter systems are theoretically simpler to 
treat and experimentally they allow much more control of different parameters. Of 
course, they cannot give any direct information about particle physics, but they 
will, nonetheless, help us understand the non-equilibrium dynamics of gauge field 
theories and in particular how they can be studied theoretically. 

Indeed, the need for reliable tools for studying non-equilibrium gauge field the- 
ories is now greater than ever before, as the Relativistic Heavy Ion Collider, which 
has already started operation, and the Large Hadron Collider, which is scheduled 
to start by 2005, will be able to study directly the quark-gluon plasma phase of 
quantum chromodynamics. The heavy-ion collisions in these experiments are so 
complicated that very reliable and accurate theoretical calculations are needed in 
order to confront the experimental results, but our present understanding of the 
theory is too rudimentary for that. The insight provided by condensed matter ex- 
periments is therefore likely to be extremely useful. In particular, it is believed 
that at a certain value of the beam energy, the quark-gluon plasma produced in the 
collision cools through a second-order transition point. ^ The relevant dynamics can 
then be described as freezing out of long-wavelength modes, ^'^ very much like in 
the theories of defect formation. 

So far, the formation of topological defects has been studied in liquid crystal^'^ 
and superfluid experiments, -"^"'-'^"'^'"'^^'^'^ which are systems with global symmetries. 
The theoretical scenario that is believed to be applicable in this case is the Kibble- 
Zurek mechanism, ^'^ but the experiments have produced mixed results. From the 
viewpoint of particle physics and cosmology, theories with global symmetries are less 
interesting than gauge field theories, whose natural condensed matter realization is 
superconductivity. The dynamics of gauge field theories are more complicated, but 
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a theory of defect formation has been developed recently-^^ for Abelian systems, such 
as superconductors. At the same time, the rapid progress in experimental techniques 
in the recent years has made accurate superconductor experiments possible. 

The purpose of this paper is to review the current theoretical status of defect 
formation in field theories, with a special emphasis on gauge field theories. The 
various applications of topological defects in other fields such as cosmology are 
outside the scope of this article, and have been discussed in other reviews. ^'^^'^^ 

Throughout this paper, we shall use natural units, i.e., c = h = Ub = ^J■o = 
The structure of this paper is the following: In Section 2, we discuss the classifica- 
tion and basic properties of topological defects in field theories. In Section 3, we 
present the basics of finite-temperature field theories and phase transitions. Defect 
formation in first-order phase transitions is discussed in Section 4. Second-order 
phase transitions are more complicated but also more relevant physically, and are 
discussed in Section 5. In Section 6, we present different approximations that have 
been used to study the dynamics of phase transitions theoretically, and in Section 7, 
we discuss different condensed matter experiments that have been carried out or 
proposed. 

2. Topological defects 

2.1. Spontaneous symmetry breakdown 

In most field theories, the Lagrangian £ [</)] is invariant under certain transformations 
(j) g{(f>), i.e., 

C[g{<P)]=C[<P], (1) 

where g is an element of the symmetry group G, typically a Lie group. The fields of 
the theory, which are here denoted by 4> for simplicity, can be in any representation 
of the group, and this determines the action of the transformation g on the field 
configuration. 

In particle physics, these symmetries are often gauge symmetries, which means 
that the group element g can be a function of space and time. For instance, the 
Lagrangian of quantum electrodynamics is 

£ = - Jf^.F'^'' + i^ii^D^, -m)i;, (2) 

where Z)^ = 9^ + ieA^^ and — d^^Ai, — di,A^ are the covariant derivative and the 
field strength tensor, respectively. This Lagrangian is invariant under U(l) gauge 
transformations 

ip{x) ^ exp{ia{x))il;{x), ^ ^^(a;) - e~^9^Q;(a:). (3) 

In the limit e ^ 0, the gauge field decouples, and the Lagrangian is only invariant 
if a is constant in space. In this case, the symmetry is said to be global. 

Gauge symmetries are extremely important in particle physics, because besides 
electromagnetism, weak and strong interactions are described by SU(2) and SU(3) 
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gauge groups, respectively. Global symmetries, on the other hand, are more common 
in condensed matter physics. 

There are situations in which a symmetry of the Lagrangian is not reflected by 
the state of the system, and the symmetry is said to be spontaneously broken. Let 
us, for instance, consider a model with a single complex scalar field 



This Lagrangian is invariant under global U(l) transformations (j) exp{ia)(p, 
where a is a constant in space-time. In a classical field theory, the vacuum state is 
simply the one in which (f> minimizes the potential V and is constant in space and 
time. If, for instance, the potential has the form 



and < 0, there is a set of degenerate minima at \(j>\ = (— m^/2A)^/^ = v, and 
therefore the system has many vacua, which arc related to each other by symmetry 
transformations. Small perturbations around any of these vacua do not possess 
the U(l) symmetry, which means that the symmetry is broken. In the presence of 
quantum or thermal fluctuations, the situation is more subtle, and we cannot say 
whether the symmetries arc broken simply by looking at the shape of the potential. 
These cases are discussed in more detail in Section 3. 

If the symmetry group G that is broken is more complicated than U(l), it may 
also be partially broken. This happens if the vacuum remains invariant under a 
subgroup H of G. Suppose, for instance, that the full symmetry group is SU(2), 
and it is broken by an adjoint scalar field $ = X]j^*cr% where a' are the Pauli 
matrices. Under symmetry transformations, $ transforms as $ — > g^^g, where 
g is an SU(2) matrix. Because all vacuum states are equivalent, we can assume 
$ = ^^a^. Then $ is invariant under transformations g = e'xp{iaa^), which form 
a U(l) subgroup of the full SU(2) symmetry group. 

2.2. Classification of defects 

Let us now consider static classical solutions when a symmetry is spontaneously 
broken. To be a solution, a field configuration must approach the vacuum asymp- 
totically at infinity. However, instead of a unique vacuum, the system has a set A4 

of vacua, which is called the vacuum manifold. In the above examples of U(l) and 
SU(2) symmetries, M is topologically a circle and a sphere, respectively. 

It is possible to have field configurations in which approaches a different point 
in A4 in different directions. In this case, the values of </> at spatial infinity form 
a mapping (f) : S^^^ M, where D is the dimensionality of the space. The 
homotopy classes of this mapping consist of configurations that can be continuously 
transformed into one another and form the {D — l)th homotopy group nn-iiM.) of 
M. Furthermore, by restricting the mapping to subsets 5" of the S^~^ at infinity, 
where n < D — 1, we obtain the lower homotopy groups 7r„(7Vt). 



c = d^cj>''d''ct,-vm. 



(4) 



vm=w?\ct>\''+x\4>\\ 



(5) 
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Because the elements of the homotopy groups are invariant under continuous 
deformations, a field configuration that corresponds to a non-trivial element of 
any 7r„(A^) with n < D — 1 cannot be continuously transformed into the vacuum 
solution. Within each homotopy class, there is a field configuration with minimal 
energy, and that configuration is necessarily a static solution of the field equations. 
However, it has typically a higher energy than the vacuum. In £) = 3 dimensions, 
configurations with non-trivial 0th, 1st and 2nd homotopies correspond to domain 
walls, vortex lines (strings) and monopolos, respectively. 

As an example, let us consider the U(l) symmetric theory (4) and a field con- 
figuration of the form (/)(r, (p, z) = \4>{r) \ exp(i(p) in cylindrical coordinates where </? 
is the azimuth angle. A configuration like this cannot be continuously transformed 
into a vacuum solution, because the path of at the infinity corresponds to a non- 
trivial element of the first homotopy group 'K\{M.). If we want (j) to be continuous, 
it must vanish at r = 0, and this costs energy. Therefore we have a classical solution 
that corresponds to a line-like object, a vortex line. 

The topology of the vacuum manifold follows directly from the symmetry break- 
ing pattern. Because different vacua are characterized by different values of M. 
is homeomorphic to the set G{<f)) of all the possible values that can be obtained 
from (j) by symmetry transformations. However, since the transformations in the 
unbroken subgroup H leave unchanged, M. is homeomorphic to the coset space 
G/H rather than the full group G. Thus, the homotopy groups TTn{G/H) tell what 
types of topological defects exist in the theory. 

In smooth configurations, which are not necessarily solutions of the classical 
field equations, we can identify the topological defects by calculating the element of 
T^ni-M) to which the behaviour of the field on a closed surface of dimensionality 
n corresponds. If this winding number is non-trivial, the surface must enclose a 
topological defect. In the above example with G =U(1), the winding number along 
a closed curve C would be 



where 9 = arg is the phase angle of the order parameter field 
2.3. Defects in gauge theories 

In gauge field theories, the discussion of the vacuum manifold becomes more subtle. 
A gauge transformation can be used to move the field anywhere in Ai at each 
point in space-time independently of all other space-time points. For instance, it is 
always possible to carry out a gauge transformation into the unitary gauge in which 
the direction of the order parameter is the same everywhere. This does not by any 
means imply that there cannot be any topological defects in gauge theories, because 
the transformation also changes the gauge field. In fact, topological defects typically 
manifest themselves in the unitary gauge as non-physical gauge field singularities 
similar to the Dirac string in the Dirac monopole configuration. 
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C = ~-F^,F^''^ + D^^*D^<jy-V{\ct>\). (7) 



A simple example of a defect in a gauge theory is the Nielsen-Olesen vortex^'^ 
in the Abelian Higgs model. This model is a gauge invariant generalization of the 
global U(l) theory in Eq. (4), and its Lagrangian is 

1 

The Nielsen-Olesen solution can be written in the form 

(f){r, if, z) = ve'^'^'ff (A/e^, evr) , A\r, ip, z) = -^tp'a{evr), (8) 

where v is the minimum of y(|0|), Nw is an integer and a{x) and /(x) are functions, 
which must be determined numerically. They vanish at a; — > and approach unity 
exponentially at a; ^ oo. If we calculate the total magnetic flux carried by the 
vortex, we find 

lim </ dr-A = - lim 2nR^a{evR) = -Nw — , (9) 
-R— >c» R—foo eH e 

which means that the magnetic flux is a multiple of the flux quantum $o = 27r/e. 
The Abelian Higgs model is very similar to the Ginzburg-Landau theory of super- 
conductivity,^® and in that context, the Nielsen-Olesen vortices are known as Abri- 
kosov flux tubes. 

An important property of the Nielsen-Olesen solution is that the energy density 
decreases exponentially at large r. In the global theory (4), the energy density 
around a vortex has a power-law decay because of the gradient term, but 
in the gauge theory it is mostly cancelled by the gauge field contribution to the 
covariant derivative Dcj). Therefore, the energy per unit length of a vortex, i.e., 
the tension, is finite in the gauge theory but logarithmically divergent in the global 
theory. This has the further implication that the interaction of two gauge vortices is 
exponentially suppressed at long distances, whereas global vortices have long-range 
interactions. 

The interaction between two vortices is attractive if A/e^ < 1/2 and repulsive 
if A/e^ > 1/2.-'^* These cases correspond to typc-I and type-II superconductors, 
respectively.^® Consequently, a solution with \Nw\ > 1 is only stable if A/e^ < 1/2. 
In the special case A/e^ = 1/2, which is known as the Bogomolnyi point, one can 
find the exact form of the functions f{x) and a{x), and calculate the vortex tension 
exactly. 

Another important defect solution is the 't Hooft-Polyakov solution in the Georgi- 
Glashow model. ^^'^^ The theory consists of a Higgs field (j) in the adjoint represen- 
tation of the SU(2) gauge group. The Lagrangian can be written as 

C = -iTVF^.F'*^ + Tr[£)^, - m^TV,^^ _ ^r,^^4^ (^q^ 

where D^^ = d^ + igA^^ and F^^i, — {ig)^^[D^, D^]. Classically, the SU(2) symmetry 
is broken into U(l) when is negative. The monopole solution has the form^^'^^ 

<P=^H{gvr), A=^[l-K{gvr)], (11) 
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where H{x) and K(x) are functions which must be determined numerically and 
which behave asymptotically as H{x) ~ x and K{x) 0. Going to the unitary 
gauge (p oc a^, we can see that the solution acts as a source for the magnetic field 
that corresponds to the residual U(l) symmetry.^^ 

Similar monopoles are also present if the group G that is broken is more compli- 
cated but has an SU(2) subgroup that is broken into U(l). In particular, practically 
any grand unified theory (GUT) predicts the existence of magnetic monopoles. 
The masscis of these monopoles would be around the GUT scale, which explains 
why they have not been observed in experiments. However, if the temperature 
of the universe was initially above the GUT scale, monopoles should have been 
formed in the phase transition in which the GUT symmetry group broke down into 
SU(3)xSU(2)xU(l). This would have had disastrous effects for the later evolution 
of the universe. Therefore, it is commonly believed that if monopoles were formed, 
they were diluted away by a subsequent period of cosmological inflation.^^ 

3. Field theories at high temperatures 
3.1. Equilibrium 



In Section 2, we discussed only classical solutions and smooth field configurations, 
but in general, we shall be interested in phase transitions that start from a thermal 
equilibrium state. The density operator in thermal equilibrium at temperature 
T=l//3isp(x exp{—(3H), where H is the Hamiltonian. From this, we can in 
principle calculate the expectation value of any observable X, 

{X) = TrpX. (12) 

For practical calculations involving only static equilibrium quantities, it is often 
easier to use the imaginary time formalism, in which the observables are expressed 
as Euclidean path integrals^^ 



-1 J £>(/.exp I- dT J d^x£Enci[Hr,x)]\ X[<^(0,f)] (13) 



and the partition function Z normalizes the result is such a way that (1) = 1. 
The path integral has a similar form in gauge field theories as well, although the 
temporal component Ao of the gauge field needs a special treatment. 

In the high- temperature limit (3—^0, the compact imaginary time dimension 
shrinks to a point. Therefore, at high temperatures, the state of the system is well 
approximated by a three-dimensional path integral^'' 



{X) = J Dcpexp (^-/3 J d-'xC3B{<Pix)]j X[0(f)]. (14) 

This approximation is called dimensional reduction, and it has been used very suc- 
cessfully in the study of the electroweak phase transition. The resulting path 
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integral has exactly the same form as the corresponding expectation value in clas- 
sical statistical physics if we replace by the classical Hamiltonian. 

In thermal equilibrium, the state of the system is specified by the expectation 
values of all the observables. In a finite volume, these expectation values possess the 
same symmetries as the Lagrangian, and therefore all observables that transform 
covariantly under symmetry transformations vanish. However, in the thermody- 
namic limit, i.e., in infinite volume, they may obtain non-zero expectation values, 
which signals that the symmetry is spontaneously broken. It depends on the param- 
eters whether this happens, and the partition function of the system is non-analytic 
at the boundary between unbroken and broken symmetries. In other words, the 
system has a phase transition. 

In the classical zero-temperature case discussed in Section 2, it was easy to 
determine whether a symmetry was broken by simply finding the global minimum 
of the potential. In the thermal case (or in the zero-temperature quantum case), 
we can define an effective potential,^° which has the same property, by 

T\n{S{^-<p)) 

VeflfM = y , (15) 

where (f) = J Sxcf). This corresponds to the free energy density of a micro- 
canonical system in which the volume average of the order parameter 4) is fixed 
to the value ip. Because the free energy is minimized in thermal equilibrium, the 
minimum of the potential shows the equilibrium state of the system. 

In perturbation theory, the effective potential can be calculated as a sum of 
one-particle irreducible (IPI) vacuum diagrams expanded around the desired value 
of (j).^^ In the global quantum theory (4), the one- loop effective potential is 



Kff('P) = + ^ 



Y,mj{^) + 0{m{^f /T) 



= V{<fi) + l\T^\^f+'^+0{m{<fifT). (16) 

This contribution increases the effective mass of ^ from to -|- AT^/3. There- 
fore, even if the symmetry is broken at zero temperature, i.e., w? < 0, it will be 
restored when the temperature becomes high enough T > Tc w (— 3m^/A)^/^. 

In gauge field theories, it is straightforward to show that the gauge invariance is 
never spontaneously broken, "^^ but this does not mean that there cannot be phase 
transitions. In a fixed gauge, we can calculate 14ff(v) perturbatively. If the Higgs 
self-coupling A is much less than the gauge coupling g"^ , this computation is believed 
to be reliable and produces a non-analytic cubic term cx jcp^.'^^''^^ This means 
that at a certain temperature T^, there are two degenerate vacua corresponding 
to zero and non-zero values of ^p. Consequently, the expectation value of (j) jumps 
discontimiously from zero to a non-zero value when the temperature decreases, i.e., 
there is a first-order phase transition. 
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When g^, perturbation theory becomes unreUable, and the hne of first-order 
transitions ends in a second-order point. In some cases, the transition disappears 
completely at higher A, and the phases are smoothly connected. In others, a 
second-order transition survives. This happens, for instance, in the Abelian Higgs 

model; the transition at large A/e^ has its most natural description in terms of a 
dual theory, whose fundamental degrees of freedom are Nielsen-Olesen vortices. ^^'^^ 

3.2. Non- equilibrium 

Formally, the time evolution of a quantum field theory is given by the operator 
equation of motion 

i^ = [X,H]. (17) 

If iJ is a constant in time, the system remains in thermal equilibrium, and if H 
depends on time, we obtain non-equilibrium evolution. The problem is that X is an 
operator, and therefore it is extremely difficult to solve Eq. (17) in any non- trivial 
field theory. We have to resort to approximations, which are discussed in more 
detail in Section 6. 



4. First-order phase transitions 

If the transition is of first order, the symmetric phase remains metastable even 
slightly below the critical temperature. The effective potential has two minima, 
and at the transition point, the one corresponding to the broken phase becomes the 
global minimum. After that, there is a finite probability per unit volume that the 
thermal fiuctuations throw the system over the potential barrier from the metastable 
symmetric phase into the true, broken minimum. It is more probable that this 
happens in only a small region of space. The boundaries around this bubble of 
broken phase cost energy, but the potential energy inside it is lower than in the 
symmetric phase. If the bubble is large enough, the contribution from the potential 
energy becomes larger than that from the bubble wall, and the bubble starts to 
grow. 

We can estimate the nucleation rate of these bubbles by calculating the energy 
Ec of a critical bubble, i.e., the lowest-energy bubble that starts to grow instead 
of shrinking away.'^^'"^*''^^ This gives the exponential suppression of the nucleation 
rate, Fnuci exp(— /3_Ec), and therefore a reasonable ordcr-of-magnitude estimate 
as long as (}Ec » 1. Calculating the prefactor requires significantly more work.^° 

In a large system, many bubbles are nucleated at different points in space, 
and when the bubbles grow, they eventually coalesce and fill the whole space, thus 
completing the transition. The bubble walls arc typically moving relatively slowly*^ 
and are preceded by shock waves, which heat up the system. This suppresses the 
further nucleation of bubbles nearby. 



4.1. Global theories 
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When a symmetry breaks, no point on the vacuum manifold Ai is preferred over 
any other. In a first-order transition, it is therefore a very reasonable assumption 
that the direction of symmetry breaking is uncorrelated between different bubbles. 
Inside a single bubble, on the other hand, it is energetically preferable for the order 
parameter field (j) to be constant. If the vacimm manifold A4 is not connected, i.e., 
T^oi-M) 7^ 1, may belong to different connected parts inside two neighbouring 
bubbles, and when the bubbles meet, a domain wall is formed between them.^^ 

If A4 is connected but not simply connected, i.e., 7Ti{Ai) ^ 1, then cj) can always 
relax to a constant after a collision of two bubbles. However, if a third bubble 
hits the first two before this has happened, a vortex can be formed. Imagine, for 
instance, that G=U(1) and inside the three bubbles the phase angle 9 of the order 
parameter (f> has the values 9i, 62 and 6^. It is common to assume the geodesic 
rule,^''^^ i.e., that when the phase angle interpolates between the values inside two 
bubbles, it always uses the shortest path in the vacuum manifold. It may then 
happen that the path that connects 9i to 92 to ^3 and finally to 9i winds around 
the vacuum manifold, in which case the order parameter cannot relax to a constant 
everywhere and a vortex is formed. Similarly, we can see that in theories with 
772 (A^) 7^ 0, collisions of four bubbles may lead to formation of monopoles. 

In this simple picture, the probability of forming a topological defect in a bubble 
collision depends essentially only on the geometry of the collision and on the geom- 
etry of the symmetry group. Therefore, we can conclude that the typical distance 
between the vortices is proportional to the typical separation between the bubble 
nucleation sites c?nuci- We can characterize the density of topological defects whose 
codimensionality is Deo by the number density n of points where a defect crosses a 
subspace of dimensionality Deo- According to the above discussion, it behaves as 

n^d-^^". (18) 

In practice, Deo is 1 for domain walls, 2 for vortices and 3 for monopoles. For 
instance, the number density of vortices per unit cross-sectional area is predicted 
to be n « d-^^y 

The probabilities of forming defects in symmetric collisions of bubbles have been 
calculated by Prokopec."'^ Although bubble nucleation is essentially a random event 
and therefore the bubble nucleation events are more or less evenly distributed in 
space and time, it was customary in early studies to assume that all bubbles were 
nucleated simultaneously in a regular lattice.*'* Obviously, the form of the lattice 
then affects the results."'^ One attempt to avoid this problem was the "lattice-free" 
approach*^'*^ in which the bubble locations are random. 

For a more accurate treatment, the field dynamics must be taken into account. 
If the bubbles do not collide simultaneously, the probability of forming a defect 
decreases, because the order parameter may already have equilibrated before the 
third bubble hits the first two. Furthermore, the growth of the bubbles is typically 
much slower than the speed of light because of the friction caused by the hot plasma 
that fills the space at high temperatures."^* This allows more time for the phase 
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equilibration between the first two bubbles to happen and therefore decreases the 
probability of forming a defect even further, although it also means that a larger 
fraction of space is in the symmetric phase at any given time and therefore the total 
number of bubble nucleation events is larger. 

4.2. Gauge field theories 

In gauge field theories, the above picture cannot be used directly, because the 
direction of the order parameter is not gauge invariant. In particular, we can always 
choose a gauge in which cj) is equal in the two colliding bubbles. That led Rudaz 
and Srivastava^^ to question the applicability of the geodesic rule in gauge theories. 

The first detailed studies of the field dynamics in the Abelian Higgs model^'^'^^ 
showed that the geodesic rule still holds if the phase transition starts from a classical 
vacuum state. In order to treat the phase angles in a gauge-invariant way. Kibble 
and Vilenkin^'^ introduced a gauge-invariant phase difference 

Agi6» = j df- DO, (19) 

where DO = V6 + eA. This cannot be used as a direct substitute for the phase 
angle difference defined in Eq. (6), because it depends on the path along which 
the difference is calculated and the phase difference around a closed curve is not 
quantized. Furthermore, Agj^ aroimd a classical vortex configuration is zero, not 
27r as in the global case, and therefore the gauge-invariant phase difference by itself 
cannot be used to locate vortices. 

Let us first consider a collision of two bubbles in a classical vacuum, where the 
magnetic field vanishes, i.e., B = V x A = 0. In this case, we can choose the 
gauge A = 0, which means that the gauge invariant phase difference is simply the 
ordinary phase difference. If we denote it by a circular magnetic flux ring with 
flux 4> = J dS ■ B = 9o/e is formed around the point of their collision, because the 
phase difference induces an electric current from one bubble to another^^'^^ (see 
Fig. la). In a simultaneous collision of three bubbles, all of the three collisions 
are at first totally independent and in each of them a flux ring like; this is formed. 
Before the collision is complete, a hole remains at the middle of the collision, and the 
total flux $tot through this hole is the sum of the gauge invariant phase differences 
divided by e, i.e., $tot = Agi^/e. When the hole shrinks, this flux is squeezed into a 
Nielsen-Olesen vortex of winding Nw = e<I>tot/27r = Agi6'/27r. Because in the ^ = 
gauge, A9 = Agi^, this is equal to the result (6) predicted by a naive application of 
the geodesic rule. 

The above discussion assumes that the fields equilibrate immediately when the 
bubbles collide, which is not always a good approximation. It was shown by 
Copeland and Saffin^"' that in certain cases the geodesic rule can be violated by 
violent non-equilibrium effects after the collision. This happens, for instance, if the 
gauge-invariant phase difference Oq is large and the bubbles are far apart, even if 
the transition starts from classical vacuum. 
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Figure 1: Bubble collisions in the Abelian Higgs model, a) In the classical vacuum 
A = 0, a, collision of two bubbles with a phase difference Oq forms a ring of magnetic 
field, b) In a collision of three bubbles, the flux rings and the initial flux $t between 
the bubbles contribute to the total flux which turns into vortices. 

However, in practice the transition starts from a thermal state, and Kibble and 
Vilenkin^^ pointed out that then the fluctuations of the magnetic field are always 
present. These fluctuations change the above picture, leading to a violation of the 
geodesic rule. More precisely, if the area between the bubbles is A, the typical 
magnetic flux will be <I>t ~ ^ AT. Suppose for simplicity that the flux is peaked 
around the centre of the collision so that it does not affect the individual two-bubble 
collisions. Then, the three flux rings are again formed, and they merge with the 
initial flux $t, which is trapped between the bubbles (see Fig. lb). This leads to 
the formation of a vortex with winding 

Nw = ( Agi6' - e$T) /27r = - e$T /27r, (20) 

where N^y is the winding number predicted by Agi^ alone. Note that Ny^^ is typically 
not an integer, but Nw always is. If the bubbles are far apart, $t can be very large, 
and therefore vortices with very high windings can form. 

The generalization of the above picture to non-Abelian gauge fleld theories is 
still an open question. It is not at all obvious what the analogue of the magnetic 
flux $T in the case of, say, 't Hooft-Polyakov monopoles is. Furthermore, because 
of confinement, a non-Abelian gauge field theory cannot be approximated by a free 
field theory. In this paper, we shall therefore only consider Abelian gauge field 
theories. 

5. Continuous phase transitions 

The dynamics of first-order phase transitions are simplified by the hierarchy between 

the bubble nucleation rate and other time scales. Furthermore, because of the 
metastability we were allowed to assume that the temperature is initially below Tc 
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and that it remains constant during the phase transition. In contimious transitions, 
the dynamics of the transition depends on how, and especially how rapidly, the 
phase transition line is crossed. 

A continuous phase transition is usually of second order, which means that the 
second derivative of the free energy is discontinuous and at least one correlation 
length diverges. Most commonly, these transitions are associated with a spon- 
taneous breakdown of a symmetry, but opposite examples are also known. For 
instance, in the two-dinicnisional XY model,'"'''' the classical analysis a, la Section 2.1 
would suggest that the symmetry gets broken, but a proper treatment shows that 
this is not the case, because a continuous symmetry cannot be broken in a two- 
dimensional theory.^^'^'' A similar theorem applies to gauge field theories in any 
number of dimensions, and therefore depending on the model, the phase transi- 
tions can range from first-order transitions^^ to smooth crossovers.^'' 

In the setup we shall consider, the system is initially in thermal equilibrium 
in the symmetric phase. The system is forced to undergo a phase transition by 
changing one of the parameters continuously through the phase transition line at a 
rate characterized by the quench timescale tq. In different cases, the parameter that 
is changed can be, for instance, a coupling constant, the mass parameter, the scale 
factor of the universe or the pressure. Although the precise way the parameters 
are changed certainly affects the detailed dynamics, the qualitative picture and the 
order-of-magnitude estimates are believed to be unaffected. 

The defects are counted after the transition when they have become well-defined 
objects. If the interaction between the defects has a short range, this can be done 
very reliably by letting the system equilibrate locally to the low temperature so 
that thermal fluctuations are negligible and counting the defects only after that. 
When the interaction has a long range, the number of defects after the transition 
has a certain ambiguity depending on how they are counted, but because we are 
only interested in orders of magnitude, we shall ignore these problems. 

In most theoretical studies, it is assumed that the phase transition is homoge- 
neous, but in a real experiment any method of cooling or changing some external 
parameters in some other way always leads to inhomogeneities, which may be im- 
portant for defect formation. ^^'^^ We shall not discuss these complications here. 

5.1. Global theories 

The overall picture of defect formation in second-order transitions in global theories 
is in many ways similar to the first-order transitions discussed in Section 4. After 
the transition, the order parameter can only be correlated at distances less than 

some finite correlation length ^ and therefore we can replace the bubbles in the 
analysis with correlated domains of radius ^.^"^'^ When we replace dnuci by ^ in 
Eq. (18), we find 

n«|-°'=°. (21) 

In cosmology, wc can state the ultimate upper limit for ^ as the age of the universe,^^ 
but in practice ^ is always much shorter than that. 
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The early suggestion^ was to identify ^ with the correlation length at the Ginz- 
burg temperature Tq defined by 



Tg 



aTcr = 1- (22) 



Below Tg, the thermal fluctuations become incapable of restoring the symmetry in 

patches larger than the correlation length. This would mean that the final number 
density of defects is independent of the quench timescale tq. However, symmetry 
restoration in a patch of radius ^ can only form vortex loops of radius less than ^ 
or monopole-antimonopole pairs with separation less than ^. and these are unstable 
configurations which quickly annihilate when the symmetry gets broken again. The 
defects that survive after the transition must have therefore been formed at an 
earlier stage. Nevertheless, the Ginzburg temperature still has the significance that 
individual defects can only be identified reliably at T < Tg- 

The attempt to use the Ginzburg temperature to estimate the defect density also 
ignores the fact that reordering the order parameter at long distances takes time.^° 
Therefore, even if it were energetically possible, it may not have have time to happen 
in a rapid transition. Instead, the relevant length scale is actually determined as a 
compromise between the tendency of ordering the fields at long distances and the 
finite time available for that process. This argument was used by Zurek to derive a 
more realistic estimate for the defect density.^'^^ 

A characteristic property of second-order phase transitions is that the correlation 
length of the order parameter diverges at the transition point with some critical 
exponent f, i.e., ^{T) « ^o|e|~'^) where e(T) = (Tc — T)/Tc and the exponent v 
depends only on the universality class of the model. However, in practice, ^ cannot 
grow arbitrarily fast, because at least it is constrained by causality to change slower 
than the speed of light. ^° If the transition takes place at a finite rate, say, e = t/rg, 
then eventually a point is reached at which the true correlation length cannot keep 
up with the equilibrium correlation length £,{T) at the same temperature. After 
that, the dynamics of the system cease to be adiabatic. 

In reality, other effects make the maximum growth rate of the correlation length 
much less than the speed of light. Zurek^ assumed that the relevant dynamics of the 
system can be characterizcxl by a rc4axation timescale t{T), which also diverges at 
the transition point in equilibrium, but typically with a different critical exponent, 
r«ro|e|-'^. 

Zurek's estimate is based on the idea that when the transition is approached 
from the symmetric phase, the dynamics of the system eventually become too slow 
to stay adiabatic and the maximum correlation length ^ is roughly the correlation 
length at that time. More precisely, he used the condition T{t) = \t\ to signal the 
breakdown of adiabaticity, which implies 

£«-(Tor^)i/(i+''). (23) 
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The maximum correlation length is then 




) 



(24) 



The critical exponents /i and u depend on the system, and at intermediate 
length scales, they may be different from their true, asymptotic values. Laguna and 
Zurek^^ have discussed two special cases, the overdamped case with /x = 1 and the 
underdamped case with ^ = 2. In mean- field theory v = 1/2, implying 



In any case, causality constrains ii>v. 

Finally, let us clarify the meaning of the correlation length by discussing it from 
another viewpoint. If we denote the dimensionality of the order parameter by N , 
there are N orthogonal elementary excitations. In the symmetric phase, they are 
all degenerate, i.e., have the same correlation length. Therefore it makes sense to 
talk about "the correlation length" ^. However, in equilibrium in the broken phase, 
only one of these correlation lengths is finite, namely the one that corresponds to 
the modulus of the order parameter. The direction of the order parameter becomes 
a Goldstone mode with an infinite correlation length. 

After a rapid phase transition, both correlation lengths are initially roughly 
equal, ^modj which corresponds to the modulus of the order parameter, and ^gs, 
which corresponds to the Goldstone modes. However, ^,„od starts to decrease 
rapidly, and because it determines the size of the defect core, the defects become 
well-defined objects. On the other hand, ^gs characterizes the separation between 
the defects, since it is easy to see that defects cause disorder to the order parameter 
field at distances larger than their separation, but at shorter distances the direc- 
tion of the order parameter is correlated. When the defect configuration evolves, 
defects annihilate, and ^gs approaches its infinite equilibrium value, albeit very 
slowly. This is another way of seeing why the maximum correlation length should 
determine the density of defects. 

In the following we shall call the above scenario the Kibble-Zurek mechanism. 
It is believed to be generally valid in phase transitions in which a global symmetry 
is broken. During the recent years, it haw been tested extensively in numerical 
simulations and condensed matter experiments, which we shall discuss in more 
detail in Sections 6 and 7, respectively. 

5.2. Gauge field theories 

As in the case of first-order phase transitions, the gauge fields were long thought 
to be irrelevant for defect formation. This view was even supported by early 
numerical simulations,^^ which showed a good agreement with the Kibble-Zurek 
mechanism. It is now understood that this was because the initial temperature was 
very low in them.^^'^'' 




(overdamped) 
(underdamped) 



(25) 
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Just like in the case of first-order transitions, thermal fluctuations of the mag- 
netic field are present at any non-zero initial temperature, and they have an im- 
portant effect on the dynamics. As discussed in Section 4.2, the magnetic flux 
that originates in thermal fluctuations can get trapped between the bubbles in 
three-bubble collisions in a first-order phase transition, leading to the formation 
of vortices with higher windings. ^'^ In a similar way, Zurek^^ argued that when a 
superconductor loop is rapidly quenched into the superconducting phase, magnetic 
flux can get trapped inside even wlicni the transition is continuous. In that case, it 
cannot be interpreted as a vortex, though, because its core is outside the system. 

The dynamics are more complicated in the "bulk" case, in which the system fills 
the whole space. It was shown by Hindmarsh and Rajantie^^ that even in a continu- 
ous phase transition, the long- wavelength thermal fluctuations of the magnetic field 
freeze out in the phase transition and form vortices. More recently, similar conclu- 
sions were also reached by Stephens et al.,^^ who studied instantaneous quenches. 
We shall call this way of forming vortices the flux trapping mechanism. 

The starting point of the fiux trapping mechanism is that in the Coulomb (nor- 
mal) phase, there are fluctuations of the magnetic field with arbitrarily long wave- 
lengths. The simplest approximation is to treat the Coulomb phase as a vacuum, 
in which case these thermal fiuctuations arc simply blackbody radiation. Provided 
that the coupling to the other fields is weak, the fluctuations are Gaussian and 
uncorrected, and the width of the Gaussian distribution is proportional to the 
temperature and independent of the momentum. 

The local magnetic fleld B is given by V x A. If we deflne its two-point correlation 
function G{k) as 

{Bi{k)B^{k')) = - ^) i2nf6 (k + k') G{k), (26) 

we have in the free-field approximation 

G{k) = T. (27) 

Of course, neither the high-temperature phase of the Abelian Higgs model nor the 
the normal phase of a superconductor is exactly a vacuum, and therefore Eq. (27) 
obtains corrections, which suppress the fluctuations. Near the transition point, 
critical fluctuations also change the correlator, and the magnetic field obtains a 
non-trivial anomalous dimension. Therefore we can generally assume that in the 
limit k ^ 0, 

G{k) oc fc", (28) 

with some anomalous dimension < 77 < 2. 

In the Higgs (superconducting) phase, long-range fluctuations are suppressed, 
and the equilibrium two-point function becomes 



(29) 
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where m-y is the inverse correlation length of the magnetic field, 




When the system goes through the phase transition, each Fourier mode of the 

magnetic field must change its amplitude from Eq. (27) into Eq. (29) in order to 
remain in equilibrium. However, it has only a finite amount of time available for 
this, and if it reacts too slowly, it falls out of equilibrium. The crucial point is that 
the modes with long wavelengths react slower and this is therefore more likely to 
happen to them. We denote by T(fc) the time scale within which a Fourier mode 
with wavenumber k can change, and in the vacuum approximation it is simply the 
inverse frequency, T{k) w . More realistically, the time scale may be determined 
by the conductivity of the material, or by plasma effects, but in any case, T{k) — > oo 
as A: — + 0, which means that there is always some critical wavelength ^ so that modes 
with longer wavelengths {i.e., k<^^'^) cannot adjust but fall out of equilibrium. 

What this means is that after the phase transition, the distribution of the mag- 
netic field at distances longer than ^ is the same as it was before the transition. 
In particular, we can calculate the typical (rms) magnetic fiux through a loop of 
radius ^ using Eq. (28). The result depends on the dimensionality D of the space, 
and assuming ?7 = 0, we find^^ 

As the modes with wavelength less than ^ were able to equilibrate, the field config- 
uration inside this loop has relaxed into its minimum energy configuration, which 
is a cluster of 

« $j-/$o « |-TV2f -^/2 (32) 

vortices, each having the same sign. Dividing this by the area of the loop, we find 
the number density of vortices formed in the transition, 

N~ 

n « ^ « _^tV2|-^/2. (33) 
^2 27r ^ ^ ' 

It is important to note that the dependence on ^ is different from the prediction (18) 

of the Kibble mechanism, and that the density increases with increasing temperature 
and increasing gauge coupling constant. In particular, Eq. (33) vanishes in the 
global limit e — > 0, and only the Kibble- Zurek mechanism remains. 

It has been argued by Stephens et al.^^ that if the phase transition is very rapid, 
the picture becomes more complicated. In the Higgs phase, the magnetic field can 
only penetrate the system if it exceeds the critical field B^i. Near the transition, 
Bel vanishes, but if the quench is instantaneous, we can assume that it jumps 
discontinuously to a non-zero value. This would mean that no vortices are formed 
in regions where the flux density is lower than B^-i. 
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Figure 2: Typical form of the real-time correlator G{t, k) oc (5(0, k)B{t, —k)) 
of the magnetic field. It is a linear combination of an exponentially decreasing 
Landau damping contribution cxp(— 71,^) and an oscillatory plasmon contribution 
exp(— 7pt) cos(a;pi). At long wavelengths, 7l <C Wp,7p, ®^ and therefore 7l de- 
termines the time r ~ 7^^^ when G{t,k) vanishes. This is the time needed for a 
perturbation of the magnetic field to equilibrate. 



In order to use Eq. (33) to predict the number density of defects, we must be 
able to calculate ^. Hindmarsh and Rajantie^^ argued that this can be found by 
using an adiabaticity condition, which balances the rate at which the equilibrium 
two-point function G{k) of a given Fourier mode must change in order to remain in 
equilibrium with the dynamical time scale of the same mode, 



dlnG(fc) 



dt 



= T-\k). (34) 



Furthermore, they argued that in the relativistic case, T(fc) is determined by the 
Landau damping rate^^ 71, [see Fig. 2] 

r-i(fc)=7L(fc)«— (35) 

Assuming that G{k) is given by Eq. (29), where w? = —5m?t/TQ, we find that 
Eq. (34) becomes 

8m? 1 4P 



from which we obtain 



(36) 



iock-^ocTQ^\ (37) 

Using Eq. (33) , we can then sec that the number density of vortices per unit cross- 
sectional area behaves in three spatial dimensions as 

n oc Tq° ^ (38) 

However, this result is not valid at very large tq.^^ When k ^ e^T, non-pcrturbativc 
effects change the form of the equilibrium correlator G{k), leading to a non- trivial 
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anomalous dimension 77. Furthermore, at low k, it is conductivity a rather than 
Landau damping that determines the rate of exponential damping 71, and this 
changes its behaviour to 7^ « fc^/cr. 

It is obvious, too, that while the equilibrium equal-time correlator G{k) is be- 
lieved to be imivcrsal, the dynamics and therefore the relaxation time scale T(fc) 
are sensitive to the details of the system. This means that the power law may be 
different for, say, superconductors, although the mechanism itself is believed to be 
valid. Finally, when applying Eq. (33) to a transition in a superconducting film, we 
have to be careful, because although the film itself is two-dimensional, the magnetic 
field extends to the third dimension as well.^^ 



5.3. Spatial distribution of topological defects 

Because neither the Kibblc-Zurek nor the flux trapping mechanism can give clear- 
cut predictions for the number density of topological defects formed in a phase 
transition unless we make certain, fairly strong assumptions about the dynamics, 
the number density itself is not a very good cjuantity for testing the mechanisms. 
Instead, we should look for signals that are truly qualitatively different in the two 
mechanisms, such as the spatial distribution of vortices after the phase transition. 
This has been emphasized by Digal et al.^^ who pointed out that the Kibble mech- 
anism leads to strong negative vortex-vortex correlations. 

A simple way of seeing that the two mechanisms predict very different spatial 
distributions is to look at the quantity Ncir),^"^ defined as the average winding 
number around a circle of radius r centred at the positive vortex. Formally, we may 
write it as the expectation value 



Nc{r) = ( n{x)^ (f dx' ■ Ve(x') \ = ln{x) ( d^x'n{x') \ , (39) 

' \x' —x\=r I \ J\x'—x\<r I 
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where n{x) is the winding number density at point x. From the latter form we 
obtain 

^^^^ = 27rr (n(O)n(r)) . (40) 
dr 

Let us first consider the Kibble-Zurek mechanism. The typical separation of 
vortices is given by ^, whereby we do not encounter any vortices of either sign at 
short distances, r <C as we increase r. However, when r ^ ^, this changes, 
because we can write Eq. (39) in the form 



Nc{r) = ^(f dx'- (n(x)V6»(f )\ , (41) 

l\x'-x\=r ^ ' 



2tt 



and by the definition of ^, the fields are uncorrelated at distance r. Therefore 
{n{x)V9{x')) = 0, and consequently, Nc{r) = 0. Thus Nc{r) must be a decreasing 
function of r, as shown in Fig. 3 by dashed lines. Using Eq. (40) this implies that 
there is a negative correlation between the vortices: Wherever there is a vortex, 
there is likely to be an antivortex nearby. 
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Figure 3: a) The behaviour of Nc{r) defined in Eq. (39) as predicted by the flux 
trapping and Kibble-Zurek mechanisms, b) The form of the vortex-vortex cor- 
relation function (n(O)n(r)) as predicted by the flux trapping and Kibble-Zurek 
mechanisms. 

The distribution produced by the flux trapping mechanism is very different. As 
Eq. (32) shows, the vortices tend to form clusters of radius ^, each containing 
vortices. Therefore, at short distances, r <^ ^, each new vortex encountered when 
r is increased typically has the same sign, and consequently dNc/dr > 0. At long 
distances, r ^, the vortex distribution is the same as the initial magnetic field 
distribution, 

(n(O)n(r)) « (B(0)B(r))equiiibrium. (42) 

In the vacuum approximation (27), the magnetic field is uncorrelated between dif- 
ferent points. Therefore (n(O)n(r)) = 0, and Eq. (40) shows that dNc/dr = 0. This 
means that at distances less than ^, Nc{r) is increasing, and near ^ it reaches a 
constant value. This is shown in Fig. 3 by solid lines. Using Eq. (40), wc can also 
say that the vortices have positive correlations at short distances: Wherever there 
is a vortex, there is likely to be another one with the same sign nearby. 

Thus, we have seen that there is a clear, qualitative difference in the vortex 
distributions formed by the two mechanisms. In practice, both mechanisms are 
believed to operate, but by measuring the form of the distribution, we can find 
out which one of them is dominant. This has already been done in numerical 
simulations, which are discussed in more detail in Section 6, but it can also be done 
fairly straightforwardly in superconductor experiments. This possibility is discussed 
in Section 7.4. 

6. Theoretical methods for studying phase transition dynamics 

Ideally, the Kibble-Zurek and flux trapping mechanisms could be easily tested nu- 
merically by solving the dynamics of the phase transition in some simple, yet non- 
trivial enough quantum field theory. However, exact analytical solutions for quan- 
tum dynamics are possible only in very simple special cases. Therefore, some ap- 
proximations are unavoidable, and many different approaches have been attempted. 
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The methods are typically variants of either the Gaussian or the classical ap- 
proximation. Examples of the first group are the linear approximation, in which in- 
teractions are neglected altogether and which corresponds to tree-level perturbation 
theory, and the Hartree approximation, which is essentially a tadpole resummation 
of the perturbativc expansion. The advantage of these approaches is that wc can, 
at least in a certain sense, keep the full quantum mechanical nature of the system. 
The price we have to pay is that the self-interaction of the field is not described 
correctly, and therefore the approximations can only be valid at early times. For the 
same reason, the Gaiissian approximations do not describe the thermal equilibrium 
state correctly near the transition and can therefore at best reproduce Eq. (24) with 
the mean-field exponents rather than with the true ones. Another weakness of the 
approach is that it is difiicult to apply to more complicated theories, such as gauge 
field theories. 

Classical approximation relies on the observation that in thermal equilibrium, 
long-wavelength modes have high occupation numbers and can therefore be as- 
sumed to behave classically. The most common realization of this idea is the 
time-dependent Ginzburg-Landau (TDGL) approximation, which consists of adding 
phenomenological noise and damping terms to the equations of motion. The main 
advantage of the classical approximation is that it describes the non-perturbative 
interactions of the long-wavelength modes correctly and, therefore, it also gives a 
correct description of the equilibrium state. However, it does not include any quan- 
tum efii'ects and it is, after all, a purely phenomenological theory, which does not 
necessarily describe correctly the dynamics of the intended fundamental theory. We 
should also be careful with the renormalization of the model, because the classical 
field theory has different ultraviolet behaviour from the original quantum theory. 

The hard-thermal-loop (HTL) approximation incorporates one-loop quantum 
corrections to the classical theory in order to avoid the weaknesses of the TDGL 
approach. It requires more resources than the classical approximation but is believed 
to approximate both the non-perturbative and the quantum mechanical aspects of 
the system correctly at reasonably long distances and time scales. 

6.1. Linear approximation 

The simplest approximation for the non-equilibrium dynamics is to approximate 
the system by a free field theory,^^ i.e., to neglect all interactions. One advantage 
of this approach is that one can study the time evolution of the full density matrix 
and therefore address issues such as decoherence and classicality.™ 

Let us consider a global 0{N) scalar field theory with the Lagrangian 

1 A 

where i, j £ {1, . . . , A^}. In the linear approximation, the equation of motion for 
each Fourier mode (j)i{k) becomes 
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If < 0, the modes with k < \J—m? axe unstable and grow exponentially. This 
approximation works in an instantaneous quench with Gaussian initial conditions 
until 

-|</).(fc)|2«-m2/A. (45) 

Because the interaction term is neglected altogether, the linear approxima- 
tion does not describe topological defects at all, but it still shows how the long- 
wavclcngth modes grow during the early stages of the time evolution. Therefore, 
if we assume that the defect distribution is related to the field spectrum in some 
particular way, it can give some information about the topological defects formed 
in the transition. A popular assumption is to identify the zeroes of the field with 
topological defects. In the special case of a Gaussian field, every observable is given 
by the two-point function G(a;), and therefore the density of zeroes can also be 
written as''^'''^ 



(no) = C4 



G"{x = 0) 



(46) 



G{x = 0) 

where the constant C depends on the dimensionality of the space and on the sym- 
metry group. 

In principle, the linear approximation can also be used for a time-dependent 

m^, but in practice the quench must be instantaneous,*'^ because othcirwisc; (0,;c6/) 
reaches —m^{t)/X immediately and the approximation breaks down. Nevertheless, 
the approach has been generalized to a slowly varying by Karra and Rivers. ^^'^^ 
In the case of a dissipationless (underdamped) field, they found ^ a Tq^, and in 

the overdamped case^^ ^ a Tq'^. Both of these results agree with Zurek's mean-field 
prediction (25). 

Non-instantaneous quenches have also been studied in a stochastic (l-l-l)-dimen- 
sional theory in the linear approximation by Lythe'^^ and Moro and hytheJ^ In the 

overdamped case,^^ the density of kinks was n oc Tq^^^ , which agrees with Zurek's 
mean- field prediction (25). In the underdamped case,*^^ they found a logarithmic 

— 1/3 

correction to Zurek's result, n (x Tq ' In rg, and they were also able to calculate 
analytically the prefactors of the power law. On the other hand, Dziarmaga''® 
has argued that the tq dependence of the underdamped case should actually be 
n (X Tq in very slow quenches. 

6.2. 2PI formalism 

One widely used way of going beyond the linear approximation is to use the two- 
particle irreducible (2PI) effective action. It is a straightforward generalization 
of the ordinary effective action to the case of a field correlator. Let us, for instance, 
consider the two-point function G{x, y) = {(j)^ {x)(j){y)) , assuming for simphcity that 
the expectation value {(j)) vanishes. We introduce an external source K{x,y), 



Z[K] = Tr 



pexp j dxdyK{x,y)(j)\x)(j){y) 



(47) 
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Then the two-point function is simply given by 

The 2PI effective action T[G] is defined as a Legendre transform 

T[G] = -ih\nZ[K] - j dxdyG{x,y)K{x,y), (49) 

where K = K[G] satisfies Eq. (48). It then follows that 

dr[G] 



dG{x,y) 



-Kix,y), (50) 



and thereby the minimum of F, given by dT[G]/dG{x, y) = 0, determines the value 
of G{x, y) in the absence of an external source. This expression is useful because 
r[G] can be written as a sum over all 2PI vacuum diagrams™ and can therefore be 
calculated fairly easily in perturbation theory. 

Although the 2PI formalism is in principle exact, in practice r[G] can only be 
calculated approximatively. In perturbation theory, the lowest non-trivial contri- 
bution comes from two-loop diagrams. This two-loop approximation is equivalent 
to the Gaussian (Hartree) approximation. ^'''^^ If we consider again the 0{N) scalar 
field theory (43), the equation of motion becomes 



d^G{k) 



k^+m^ + {2 + N)xl ^G{k') 



G{k), (51) 



where we have defined G{k) as the Fourier transform of G{x) given by {(pi {x)<pj (y)) = 
dijG{x — y). As in the linear approximation, we can then use Eq. (46) to estimate 
the number density of topological defects. 

Because gauge fields lead to technical difficulties within this approach, it has 
only been applied to phase transitions in theories with global symmetries. The 
density of topological defects formed in an instantaneous quench was calculated 
by Antunes and Bettencourt,^^ who also compared their results with the linear 
theory. Later, the Hartree approximation was also applied to non-instantaneous 
quenches in a two-dimensional theory with a global 0(2) symmetry by Bowick and 
Momen,^^ but they did not calculate the density of defects explicitly. In the case of 
an expanding universe, a proper numerical test for Zurek's predictions was carried 
out by Stephens et al.^^ who found the scaling laws ^ ex Tq^^ for underdamped 
and ^ cx Tq^^ for the overdamped cases, in a reasonable agreement with Zurek's 
predictions (25). 

Although the Hartree approximation is sometimes called non-perturbative in 
the literature, it is only non-perturbative in the same sense as the calculation of the 
effective potential (15) using one-particle irreducible diagrams. This means that 
the critical exponents ji and v in Eq. (24) produced by this approximation arc not 
the true critical exponents of the system, but the perturbative ones. Therefore, the 
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Haxtree approximation, or indeed the 2PI formalism at any finite loop order, cannot 
give a correct description of the true dynamics of the phase transition. Nevertheless, 
it is still a non-trivial test for Zurek's predictions to check whether the results of 
the Hartree approximation agree with the predicted mean- field behaviour (25). 

Being perturbative, the Hartree approximation is also insensitive to truly non- 
perturbative objects such as topological defects and therefore it cannot describe the 
formation of topological defects properly. This is obvious from Eq. (51), because 
the choice of the symmetry group 0(A'') only affects the strength of the effective 
coupling between different modes, whereas in reality the dynamics at different N 
are very different: N = 1 gives rise to domain walls, A'^ = 2 to vortices and N = 3 
to monopoles. After all, the Hartree approximation is more or less equivalent to the 
leading term in the large- -/V approximation, and there are no topological defects in 
the largo- limit. Another flaw in this approximation is that it does not describe 
scattering of different Fourier modes properly, because each mode only interacts 
with the average of all the other modes. 

The reason why Eq. (51) has such a simple form is that all the fundamental 
fields in the theory are related by a symmetry and therefore they all have the same 
two-point function. In more realistic theories, the equation of motion becomes more 
complicated, and the application of this formalism to gauge field theories is largely 
an unsolved problem. A further problem is that especially in gauge theories, even 
the equilibrium properties become non-perturbative near the transition and are 
therefore not correctly described by the perturbative 2PI effective action. 

At two-loop level, there is no difference between classical and quantum mechan- 
ical time evolution. Because of this and the other problems mentioned above, it 
is essential to be able to go beyond the Gaussian approximation. At three-loop 
level, scattering of modes is taken into account, although it must still be as- 
sumed to be weak. In this approximation, some differences appear between clas- 
sical and quantum dynamics, but the equations of motion still look very similar. 
One possible alternative to the perturbative calculation of the effective action r(G) 
is the large-A'' expansion,^^ but as mentioned above, it cannot describe topologi- 
cal defects. Other ways of improving the Gaussian approximation have also been 
proposed. ^^'^"'^^ 

6.3. Classical approximation 
6.3.1. Basics 

If we are willing to sacrifice the quantum mechanical nature of the system, it is 
relatively straightforward to study the dynamics in a fully non-perturbative way 
by simply solving numerically the classical equations of motion. In fact, this is 
believed to be a good approximation as long as the temperature is high enough, 
because the defect density and other physically relevant properties of the system 
after the transition are determined by the long-wavelength field modes, and they 
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have high occupation numbers in thermal equiHbrium.^^ 

A significant advantage of the classical approach over the 2PI formalism is 
that in the case of static equilibrium observables it is equivalent to dimensional 
reduction,^^'^^ i.e., to an approach in which the full four-dimensional path inte- 
gral is replaced by an effective three-dimensional one [sec Eq. (14)]. Because the 
construction is free from infrared problems, it is believed that the effective theory 
describes the phase transition correctly. In particular, it belongs to the same uni- 
versality class as the full quantum theory, and therefore the critical exponent u in 
Eq. (24) has its true value. However, as discussed in Section 3, the imiversality 
arguments do not apply to the time dependence, and therefore the other critical 
exponent ^ is not necessarily reproduced correctly. Nonetheless, the classical sim- 
ulations provide a highly non-trivial, non-perturbative test for the mechanisms of 
defect formation. 

The general strategy in classical simulations is to take an initial field configura- 
tion that is in thermal equilibrium and follow its time evolution by solving numer- 
ically the classical equations of motion. The thermal initial condition is typically 
prepared with a Monte Carlo algorithm of some type.^'^'^'' Perhaps the most popular 
choice is to use Langevin dynamics, ^^'^^ which means adding noise and damping 
terms to the equations of motion. In the case of an 0{N) symmetric scalar field 
theory, this means 

dofpi + 'rido(t>i - ^'^<f>i - m^fpi - X(l>j(l>j(l>i = Ci, (52) 
where the white noise term ( satisfies 

{Q{t,x)Qit',x')) = 2TrjS,jS{t - t')S'^''\x ~ x') . (53) 

In the limit of infinite time, t oo. and zero timc^ stcip, the probability distribution 
of the field configurations approaches the canonical ensemble 

p[cj)\ a exp(-i?[0]/r), (54) 

where H\4)\ is the Hamiltonian. 

Another popular approach is the Metropolis algorithm, '^^ which consists of mak- 
ing small random changes in the field configuration {0-} and calculating the 
corresponding changes A£ in the total energy. If A£ < 0, the change is accepted, 
and if AE > 0, it is accepted with probability exp{—AE/T). When the number 
of these updates is large, the probability distribution of the field configurations 
again approaches the canonical ensemble. The Metropolis approach is generally 
more flexible than the Langevin approach, because the nature of the small random 
changes can be chosen to optimize the thermalization rate. In certain special cases, 
modifications of the Metropolis algorithm such as the heat bath algorithm,^® in 
which one degree of freedom is brought in contact with a heat bath at a time, and 
hybrid Monte Carlo, in which one uses the Hamiltonian time evolution to move 
the system around in the phase space, are more efficient. 
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In three dimensions, the energy density of a classical field theory in equilibrium 
is ultraviolet divergent, but this does not constitute a severe problem, because the 
space and time must be discretized in any case in order for the numerical solution 
to be possible, and this lattice cutoff A = 1/a provides the necessary ultraviolet 

rcgularization. Nevertheless, these extra contributions should be rcnormalized, i.e., 
the bare couplings used in the simulations must be chosen in such a way that they 
correspond to the desired infrared physics. In the limit of very fine lattices, o — > 0, 
we can calculate perturbatively the counterterms that are sufficient for rcnornializ- 
ing all static equilibrium obscrvablcs.^^"'^*'^ but in gauge field theories, they leave 
some lattice spacing dependence in the time evolution. In particular, this dis- 
crepancy becomes worse in the continuum limit, and therefore we should not use 
too fine lattices. 

It is quite common to keep phenomenological noise and damping terms in the 
equations of motion even during the time evolution of the system, in order to approx- 
imate the coupling to the other degrees of freedom. This time-dependent Ginzburg- 
Landau (TDGL) approach leads to a faster thermalization than the Hamiltonian 
equations of motion, but it does not change the nature of the dynamics qualitatively, 
because the interactions lead to thermalization in any case. 

One advantage of classical simulations is that we can locate the defects directly 
in the field configurations and do not have to rely on the Halperin-Mazenko-Liu 
approximation (46). In a global U(l) theory, the lattice winding number is a 
straightforward analogue of the continuum one (6). For each link {x,x + i), we 
define 

^0{x + i) -0{x)\ , (55) 



where the notation [. . .J^r indicates that the difference is calculated in such a way 
that it always lies between — tt and n. This is equivalent to the geodesic rule. The 
winding number of a plaquette is then simply [cf. Eq. (6)] 

^^^^^ = h ("^^(^.^+2) + - ^^(^+3,5+2+3) - • (56) 



6.3.2. Global simulations 

Antunes and Bettencourt^^ compared the dynamics of the classical (l-l-l)-dimen- 

sional scalar field theory in an instantaneous quench with linear and Hartrec ap- 
proximations. Their results show that both approximations break down when the 
back-reaction becomes important, i.e., when the order parameter reaches the min- 
imum of the potential. Laguna and Zurek^*^^'^"^ studied the formation of kinks in 
a linear quench in the same model in the overdamped case. They found a power 
law dependence n oc oc -7--0-28±o.02^ which agrees with the mean-field Zurek 
prediction (25). Later, they investigated different cases that interpolate between 
the overdamped and undcrdamped dynamics and found power-law exponents rang- 
ing from —0.23 ± 0.01 (overdamped) to —0.33 ± 0.01 (underdamped) . Moro and 
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Lythe*^^ also simulated the same system and confirmed their analytical prediction 

for the prcfactor of the power law. 

The reason why the above results are consistent with the mean-field expectations 
is that there is actually no phase transition at all in any (l+l)-dimensional classical 
model, and therefore no critical behaviour either. This also means that we cannot 
unambiguously measure the final defect density. In order to test the scenarios of 
defect formation properly, we have to go to higher dimensions. 

Antunes et al}^^ studied the formation of vortices in a (3+l)-dimensional scalar 
field theory with a global 0(2) symmetry. They used the TDGL approach and, 
instead of varying the mass parameter m^, they changed the amplitude of the noise 
term, which effectively means changing the temperature of the heat bath to which 
the system is coupled. The final density of vortices per unit area had a power-law 
dependence on the quench rate n cx oc 7--0-4982±o.079 ^pj^jg compatible with 
the overdamped mean- field prediction (25), although the classical simulation should 
experience the true critical behaviour with non-mean-field critical exponents. This 
is an indication that the freeze-out of the order parameter field took place before 
the system entered the critical region. 

In their simulations within the same model, Bettencourt et al.^'^'^ first carried 
out a linear temperature quench but then heated the system slightly above the 
critical temperature for a short time. They found that the original vortex network 
survives this period of symmetry restoration, which supports Zurek's idea that the 
defect distribution is determined by degrees of freedom whose dynamics are very 
slow near the transition even in the symmetric phase. 

The model was also studied by Bowick et al}'^^ who used a time-dependent 
mass term to start the phase transition. They concentrated on the evolution of the 
correlation lengths and showed that a naive application of Eq. (21) overestimates 
the number of defects by a factor of 2 ... 4, because the correlation lengths keep on 
growing after the transition. 

Stephens^°^ used the TDGL approach to study the phase transition in the (2-1-1)- 
dimensional scalar field theory with a global 0(3) symmetry. The topological defects 
in this model are textures, whose energy is independent of their size in the classical 
vacuum. Therefore not only the density but also the typical size of these objects 
is determined by the phase transition dynamics. The results show that the typical 
distance of the textures behaves as Lgep oc f^^^^^-^"^ and the size as oc ^-^■^^±0.04^ 

whereas Zurek's prediction (25) for the freeze-out scale would have been ^ oc Tq^^ . 
However, these length scales evolve dynamically after the transition, which may 
explain the discrepancy. 

6.3.3. Gauge field theories 

One clear advantage of the classical approximation is that it can be easily applied 

to models with gauge symmetries. We shall discuss only the Abelian Higgs model 
here, but in principle similar techniques can and have been applied to non-Abelian 
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theories. The classical equation of motion for the gauge field is 

D^F^''=r, (57) 

where = —2elm<p* D'^ <p is the electric current. It is most convenient to use the 
temporal gauge Aq =0, because that simplifies the equation of motion into 

doE = VxVxA + j, (58) 

where E = —doA is the electric field. However, in the temporal gauge we also 
obtain an extra constraint, which the field configurations must satisfy and which is 
nothing but the Gauss law 

V-E = p, (59) 

where p is the electric charge. 

The necessary discretization of this model is straightforward. There are two 
alternative formulations, compact and non-compact, which correspond to gauge 
groups U(l) and M, respectively. Because vortices are not absolutely stable in the 
compact formulation, we only consider the non-compact formulation here. The 
gauge field is described by a real- valued field = ae^^, which is defined on links 
between lattice sites. 

It is obvious that Eqs. (55) and (56) cannot be used in gauge field theories, 
since they are not gauge invariant. If the gauge is not fixed, the phase angles 
are totally random and uncorrelated, and therefore it is always trivially true that 
Nij = 1/3.^^^ The solution is to "relax the geodesic rule" by replacing the phase 
difference in Eq. (55) with its gauge-covariant counterpart^^^'^^^ [cf. Eq. (19)] 

- aiix). (60) 

When this is used in Eq. (56), the winding number becomes gauge invariant, and it 
is easy to see that it is always an integer. In fact, it is the natural lattice analogue 
of the continuum winding number Nw defined in Eq. (20). 

6.3.4. Gauge simulations 

Yates and Zurek^^ investigated the phase transition in the (2-|-l)-dimensional Abe- 
lian Higgs model using the TDGL approach. They measured the number density 
of vortices after the transition, and found power laws n oc oc 7-^"-44±o.io 

the overdamped case and n oc oc Tq°-^°^°-°'' for the underdamped case. These 
results agree with the mean-field prediction (25) based on the Kibble-Zurek mech- 
anism. The values they used for the gauge coupling constant and temperature, 
e = 0.5 and T ~ 0.01, were so low that the flux trapping mechanism (33) does not 
have significant effects. Indeed, the typical flux through their whole 512^ lattice 
is merely $ « 4$o according to Eq. (31). They also used Eq. (55) to define the 
winding number, but again because of the low temperature, this did not affect the 
results. 
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The (2+l)-dimensional Abelian Higgs model was also studied by Ibaceta and 
Calzetta/^'' who used the overdamped TDGL approach. Because they used an in- 
stantaneous quench, they were not able to test the predictions for the defect density. 
On the other hand, their aim was to test the applicability of the linear approxima- 
tion, and they found that although it eventually breaks down, it agrees very well 
at early times. This is presumably partly because they used Gaussian initial con- 
ditions for the Higgs field rather than thermal ones. In their initial conditions, the 
gauge field vanished everywhere, so they could not have observed the flux trapping 
mechanism in any case. Furthermore, they used Eq. (55) as the definition for the 
winding number. 

The only (3-|-l)-dimensional simulations in a gauge field theory so far have been 
carried out by Hindmarsh and Rajantic.-'^^'^^ In Rcf. 14, they used the classical ap- 
proximation without damping and noise terms. In their lattice, one of the dimen- 
sions was shorter than the other two, in order to stabilize vortices that wind around 
that dimension so that they can be counted. They prepared the initial configuration 
using a hybrid Monte Carlo algorithm, and used a time- varying mass term in their 
simulations. They used the gauge-invariant definition (60) for the winding number 
to locate the defects in the final state, and found the power-law n oc t--o.250±o.oi3 
when the short dimension was 5 in lattice units, and n a t--0-274±o.039 ^jjgj^ ^g^g 

20. These results disagree strongly with Zurck's prediction (25), which is n oc Tq 
On the other hand, they can be explained in terms of the flux trapping mechanism. 

Prom the viewpoint of vortex dynamics, the system they studied was eff'ectively 
two-dimensional and therefore the relevant theoretical predictions correspond to 
£> = 2 in Eq. (33), i.e.,^^ 

e TV2 

where is the extent of the system in the short dimension. This predicted depen- 
dence on Lz was confirmed in the simulations."'^^ Note, however, that using Eqs. (37) 
and (61), we expect n oc Tq'^'^. This discrepancy is explained by the failure of 
the classical approximation to describe the time evolution of the quantum theory 
correctly. In particular. Landau damping (35) is very sensitive to the quantum 
mechanical ultraviolet modes, and in the classical theory, it actually behaves more 
like***^ 7L(fc) oc fc^-i. Substituting this into Eqs. (34) and (61), we find n oc Tq°-^^, 
in a fair agreement with the measurements. 

As discussed in Section 5.3, the flux trapping mechanism also predicts a charac- 
teristic spatial distribution for the vortices. Hindmarsh and Rajantie^'* studied that 
in their simulations by measuring the function Nc{r) [see Eq. (39)]. Their results 
for two different quench rates are shown in Fig. 4, and they display a clear positive 
correlation at short distances, as predicted by the flux trapping mechanism. 

More recently, Stephens et al.^'^ studied instantaneous phase transitions in the 
(2-1-1 )-dimensional Abelian Higgs model using overdamped TDGL simulations. They 
foimd very strong clustering of vortices, as shown in Fig. 5, in agreement with the 
discussion in Section 5.3. They also measured the dependence of the vortex density 
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Figure 4: The correlation function Nc{r) [see Eq. (39)] measured by Hindmarsh and 
Rajantie^^ in two quenches with different tq. The solid line illustrates the effect 
of the finite system size by showing Nc{r) for a random distribution of vortices 
and antivortices. In both runs, Nc{r) is above this curve, which supports the flux 
trapping scenario, as discussed in Section 5.3 (see also Fig. 3). (Prom Ref. 14.) 




Figure 5: An example of the spatial vortex distribution after a phase transition in 
the simulations of Stephens et al.^^ White and black circles correspond to vortices 
and antivortices, respectively. The plot shows clearly the clustering of vortices 
discussed in Section 5.3. (From Ref. 64, kindly supplied by G. Stephens.) 
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on the gauge coupling and temperature, and found the relation n w CeT^/"^, where 
C is a constant. This confirms the prediction (33) of the flux trapping scenario. 

People have also studied inhomogeneous setups with the motivation that they 
might be more relevant for actual experimental phase transitions.^^ Dziarmaga et 
al}^^ used the TDGL approach in a global theory in (1+1) dimensions to study 
a phase transition which starts from one end of the system and propagates to the 
other. They found significant suppression of defect density. Aranson et al}^^'^^"^ 
have used overdamped TDGL simulations to study a global (3+l)-dimcnsional U(l) 
theory in an inhomogeneous case, in which the system was initial heated up locally 
at a single point. 

6.4. Hard-thermal-loop improvement 

When discussing the classical approximation in Section 6.3, we implicitly assumed 
that the classical theory had the same Lagrangian as the full quantum theory it was 
supposed to approximate. However, Bodeker et a/.^°^'^°^ showed that this leads to 
wrong dynamics, although it reproduces the static equilibrium properties correctly, 
once simple renormalization countcrtcrms are included. 

The reason for this failure is that the classical approximation is only valid for 
long-wavelength modes. Therefore, the classical theory must be interpreted as an 
effective theory in the Wilsonian sense, and in principle, its Lagrangian should 
contain all the terms that are compatible with the symmetries of the system. In 
vacuum, these terms are strongly constrained by the Lorentz invariance, but at a 
non-zero temperature, Lorentz invariance is broken by the rest frame of the thermal 
background. Therefore many more terms are allowed in the effective Lagrangian 
and restricting it to its original form is not justified. 

Fortunately, the difference between the quantum and classical theories appears 
only at very high momenta, and these ultraviolet degrees of freedom are pertur- 
bative. Therefore wc can calculate the necessary corrections to the classical La- 
grangian perturbativcly. In a global scalar theory, we only obtain a simple mass 
counterterm,^^^ but in gauge field theories, the correction is a complicated non-local 

object. In the Abelian Higgs model, this hard-thermal-loop improved effective La- 

• 1 2n 
grangian is 

£htl - ~F,,F - -m^ j —F j^F,^ 

+ |D^</.|2-m2|</.|2_A|</.|4, (62) 

where rn^ = rn? + (e^/4 + A/3)r^, and the integration is taken over the unit 
sphere of velocities v — (l,w), 'tP' — 1. The Debyc screening mass has the value 
m\) = |e^T^. We have here ignored a similar, non-physical, contribution that 
arises from ultraviolet lattice modes, and which should in principle be subtracted. 
Therefore wc must assume that the lattice spacing is long enough, a ^ T~^, so 
that this contribution is negligible. Because the perturbative calculation that gave 
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Eq. (62) involved only hard, ultraviolet modes, it is believed to be reliable as long 
as the coupling constants are small. 

In principle, we can simulate the dynamics of the theory by solving numerically 
the equations of motion derived from Eq. (62). However, because the extra term is 
non-Iocal this requires keeping all the previous time steps in the memory. Instead, 
it is much more convenient to introduce extra degrees of freedom, which have local 
interactions, and which reproduce the same effective term. There are different ways 
of doing this: The extra degrees of freedom can either be point particles^^^ or 

fields.124.122,125 

The HTL approach has been used to study equilibrium quantities such as the rate 
of baryon number violation in the symmetric phase of the electroweak theory.^^^'^^^ 
Similarly, Hindmarsh and Rajantic*'^ used it in the Abelian Higgs model to measure 
the real-time correlator of the magnetic field at the transition point. They found 
that within the range of wavelengths they studied, the time scale that dominates 
the dynamics of the long-wavelength magnetic fields is determined by the Landau 
damping rate 7^ and that at intermediate wavelengths it agrees with the pertur- 
bative result (35). Because 71, is very sensitive to quantum effects — they also 
measured the rate without the HTL correction and obtained a very different re- 
sult — this is a very strong indication that the HTL approach really approximates 
the full quantum dynamics well. Nonetheless, because the approximation ignores 
the scattering of hard modes, it eventually breaks down at very long wavelengths, 
fc < e^T, and very long times. 

The HTL approximation can also be used in non-equilibrium settings provided 
that the ultraviolet modes are still in equilibrium, which is the case in thermal 
phase transitions. Hindmarsh and Rajantie^^ used this approach to study vor- 
tex formation in the phase transition of the Abelian Higgs model, and found n oc 
T--0-'2oi±a.oi5 ^ This differs from the prediction of the Kibble-Zurek mechanism (25) 
which is n (X ^"■^ oc j--0-5 - -0-66^ agrees very well with the flux trapping predic- 
tions (37) and (61). 

7. Experiments 

The similarity of defect formation in condensed matter systems and cosmology was 
already noted by Zeldovich et al.^^ and Kibble, but the first concrete proposal 
of utilizing this correspondence to test the cosmological scenarios in experiments 

was made by Zurck."'''^^ He suggested an experiment in which a pressure quench 
would be used to cause a phase transition in ''He and which was later carried out 
by Hendry et al}^ 

Later on, other condensed matter systems have been used to study defect forma- 
tion as well. Most of them, liquid crystals, '^Hc^"'^'' and ■^He-B^^'^^ are systems 
with global symmetries, and therefore the applicable theoretical scenario is the 
Kibble-Zurek scenario. The same applies to recent studies of non-linear optical 
systcms^^^ and convection in fluids, as well as to the proposed experiments with 
atomic Bose-Einstein condensates. 
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The only exception are superconductors.^^ The order parameter, the Cooper 
pair, is electrically charged, and therefore the symmetry that is broken is a local 
symmetry. This means that the flux trapping mechanism should form vortices in 
addition to the Kibble-Zurek mechanism. So far, only two superconductor exper- 
iments have been carried out,^^^'^'^" and they haven't been able to give conclusive 
results, but the rapid progress in experimental techniques in the recent years sug- 
gests that more detailed experiments will soon be possible. 

7.1. Liquid crystals 

Liquid crystals^'^^ arc perhaps the simplest condensed matter systems in which 
formation of topological defects can be studied. The phase transition between 
isotropic and nematic phases takes place near the room temperature, and the defects 
can be seen with an optical microscope. A disadvantage of liquid crystals is that 
the transition is of first order and therefore it cannot really probe the Kibble-Zurek 
mechanism, which applies to continuous phase transitions. 

A nematic liquid crystal consists of rod-like molecules, whose orientation acts 
as the order parameter. It is an unorientcd three-vector and is commonly denoted 
by n. The vacuum manifold is S"^ jlii and has non-trivial first, second and third 
homotopy groups. Therefore, the system has three types of topological defects: 
vortex lines (strings), monopoles and textures. 

The first experimental study of defect formation was carried out by Chuang et 
al.^ and concentrated mainly on the dynamics of the string network after the tran- 
sition from the isotropic to the nematic phase. Bowick et al.^ paid more attention 
to the formation of the defects, and carried out a quantitative comparison with the 
Kibble mechanism. Because the phase transition is of first order, the number den- 
sity of vortex lines in the nematic phase should be given by Eq. (18). They counted 
the number of bubbles, compared it with the number of vortices and found a good 
agreement with the theoretical prediction. 

As discussed in Section 5.3, the spatial distribution of defects is a very convenient 
way to test the mechanisms for defect formation. Because the symmetry broken in 
the isotropic-nematic phase transition is global, the relevant theoretical scenario is 
the Kibble mechanism, and it predicts strong negative correlations between defects. 
These correlations were studied experimentally by Digal et al.^^ They measured 
the net number of vortices Nw{A) through a loop of area A. On average, this is 
obviously zero, because vortices of either sign have the same probability, but the 
width a of the distribution depends on the correlations of the vortices. Digal et al. 
assumed a power-law behaviour 

a = CN", (63) 

where C and u are free parameters. They used the total number of vortices N inside 
the loop instead of A, because it separates the effects of the vortex correlations from 
the relation between N and A. 
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The Kibble mechanism predicts v ~ 1/4, but if the vortices were randomly 
distributed, the exponent would be = 1/2. The measured value for the exponent 
was v = 0.26 ±0.11, in a good agreement with the Kibble mechanism. Moreover, 
Digal et al. measured the value of the prefactor C = 0.76 ± 0.21, which also agrees 
well with the theoretical prediction C w 0.71, calculated using the simplifying 
assumption that the bubble nucleation sites form a square lattice. 

7.2. Superfluids 

At low enough temperatures, helium becomes a superfluid,^^ and this transition 
can also be described as a spontaneous breakdown of a global symmetry. It has the 
advantage over liquid crystals that superfluidity is a quantum phenomenon, and 
therefore the experiments can really probe the behaviour of a quantum field theory 
rather than a classical one. 

In the '^He isotope, the phase transition to the superfluid phase takes place at 
the critical temperature T\ 2.18 K. In the Ginzburg-Pitaevskii^^^ picture of the 
transition, the order parameter ijj is the quantum mechanical wave function of the 
^Hc atoms. In practice, -0 is a complex scalar field, and the theory is invariant under 
global U(l) transformations. In the superfluid phase, the ^He atoms form a Bose 
condensate, which is signalled by a non-zero vacuum expectation value of ■0. This 
breaks the U(l) symmetry, and as discussed in Section 2, leads to the existence of 
vortex line solutions. 

Vortex formation in a ^He phase transition was first studied experimentally by 
Hendry et al}° They had a small volume of ^He in a container at a tempera- 
ture slightly above Tx, and instead of cooling the system, they expanded the con- 
tainer rapidly so that the pressure decreased. This increases Tx, and therefore 
the system underwent a phase transition into the superfluid phase. They counted 
the vortices formed in the transition by measuring the; attenuation of the second 
sound in the liquid after the transition. The results show that vortices were indeed 
formed in the transition, but this was later attributed to hydrodynamic effects that 
arose from non-idealities in the experimental setup rather than to the Kibble-Zurek 
mechanism. 

A more careful study was carried out later by Dodd et al}^'^^^ who eliminated 
the most significant sources for hydrodynamic vortex formation. Surprisingly, they 
did not find any evidence for vortex formation in the transition. They suggested 
that the vortices may have decayed faster than expected. This idea was made more 
precise by Karra and Rivers,^^ who argued that the Ginzburg temperature Tq [see 
Eq. (22)], above which thermal fluctuations are still able to restore the symmetry 
within a correlated region, is well below Tx in ^He, and in fact, the final state was 
still above Tq in the experiment. Therefore, the vortices never became well defined 
classical objects and were washed out by thermal fiuctuations before being observed. 

The other helium isotope ^He also displays superfluidity at low temperatures,^^ 
but many details are very different. Because ^Hc atoms are fcrmions, they can 
only condensate if they form bosonic Cooper pairs. Therefore, superfluidity needs 



Formation of topological defects in gauge field theories 35 



millikelvin temperatures, much lower than in ^He. Furthermore, the Cooper pairs 
may have a non-trivial spin S and angular momentum L, which means that the 
order parameter is more complicated than simply a complex scalar field. Therefore, 
there are different superfluid phases: the A phase, which is only present in a narrow 

range of temperatures and under high pressures or if a strong external magnetic 
field is applied, and the B phase. In an external magnetic field, a third superfluid 
phase, known as the Al phase, is also possible. 

In all the existing experiments, the transitions have been between the normal 
and the B phase, and therefore we shall only concentrate on that. In the B phase, 
the total angular momentum J = L + S vanishes. The symmetry breaking structure 
isi9 

S0(3)l X S0(3)s X U(l)^ - S0(3)l+s- (64) 

This allows two different types of vortices, oriented mass vortices and unoriented 
spin vortices}^^ The latter ones are bound to soliton sheets and are therefore un- 
stable. 

Vortex formation in ^He has been studied in two experiments. ^^'^^ They both 
used external neutron sources to heat up the system locally. When a neutron hits 
a ^He nucleus, it may be captured in the reaction 

n+^He^^H + p. (65) 

This releases 764 keV of energy, which heats up a "hot spot" of radius « 30 /im, 
and inside this hot spot the symmetry is locally restored. The hot spot cools 
rapidly in about 1 jus, and undergoes a phase transition back to the superfluid 
phase. According to the Kibble-Zurek scenario, a tangle of vortices is formed in 
this process. 

In the Helsinki experiment by Ruutu et al}^ the cryostat that contained the 
helium was rotating with a velocity lower than the critical velocity of the superfluid. 

This rotation expands those vortex loops that are larger than a certain critical size. 
Eventually, they straighten out and move to the centre of the container. Using 
NMR spectroscopy, Ruutu et al. were able to count the vortices one by one and 
thereby obtained a very precise measurement for the number of vortices formed in 
the transition. 

In the Grenoble experiment carried out by Bauer le el al}^ the number of vortices 
formed in the transition was inferred using a vibrating superconducting wire, which 
measures the amount of heat deposited in the liquid very accurately. If vortices are 
formed, this amount is less than the total 764 keV released in the neutron capture 
event, because some of the energy is stored in the vortices. The experiment was 
carried out at a temperature well below so that the lifetime of the vortices was 
much longer than the time needed for the measurement. 

In both "^He experiments, the measured number of vortices was in a reasonable 
agreement with the prediction of the Kibble-Zurek mechanism. In experiments of 
this type, the quench timescale tq is not a free parameter, because it is determined 
by the equilibration rate of the system after a neutron capture event, and therefore 
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it is not possible to study the dependence of the vortex number on tq. Nevertheless 
the results show beyond doubt that vortices are formed in the transition, in contrast 
to ^He, which may be explained by the fact that in ^He, Tq is very close to TcJ^ 

In the above experiments, the observed vortices were mass vortices, but more 
recently, Eltsov et al,}"^^ have also reported an observation of composite spin-mass 
vortices after a phase transition to the B phase. 

Another type of a superfluid experiment was also suggested by Zurek.^ He 
pointed out that if the superfluid container is an annulus instead of a cylinder, 
the "vortices" that arc formed in the transition pass through the hole. Because 
their cores are outside the superfluid, they cost very little energy and are therefore 
very stable. They can be detected as a non-zero angular velocity of the superfluid 
around the annulus. However, this experiment has not been carried out yet. 

7.3. Other global systems 

The liquid crystal and superfluid experiments were not able to measure the depen- 
dence of the vortex number on tq and therefore they could not test the Kibble- Zurek 
mechanism quantitatively. However, there have recently been reports of other stud- 
ies in which this has been done.^^®'^^^ 

The first such experiment was carried out by Ducci et al}^^ in a non- linear 
optical system. They studied a liquid crystal light valve illuminated by a laser 
beam and inserted in a feedback loop. When the intensity of the incident light 
exceeds a threshold value, an intensity patterns of standing rolls forms. The authors 
changed the intensity of the light at a finite rate characterized by tq , recorded the 
intensity pattern and counted the defects. They found a power-law dependence 
n oc 7--0 '''"±o.o4^ ^ very good agreement with Zurek's prediction (25). 

Casado et al}^"^ studied the breaking of a global symmetry in the Benard- 
Marangoni conduction-convection transition. They heated up the bottom of a cylin- 
drical container with a layer of silicone oil in it. Below a certain critical temperature 
Tc, the fluid is in a homogeneous conduction state, but above Tc, the symmetry is 
broken by a hexagonal array of convection lines where hot fluid flows upwards. 
Casado et al. captured the image with a CCD camera and located the defects in 
the array structure. They used different rates of increasing the temperature at the 
bottom of the container, characterized by the quench time scale tq. The number 
of defects had a power-law dependence on tq, with exponent ranging from —0.45 
to —0.25 depending on the viscosity of the fluid. 

7.4. Superconductors 

All the experiments discussed above are systems with broken global symmetries, 
whereas local gauge symmetries are more relevant for particle physics and cosmol- 
ogy. Therefore, it has been appreciated for a long time^^'®^'^^^ that defect formation 
should also be stiidicd experimentally in type-H superconductors. 

In superconductors,^^ the order parameter is a scalar fleld tp, which describes the 
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Cooper pairs. The field is charged under the U(l) gauge group of electrodynamics, 
because the Cooper pairs have an electric charge of — 2e. In the superconducting 
phase, the Cooper pairs condense, and this leads to the Meissner effect. The equi- 
librium properties of superconductors near the phase transition are described by 
the Ginzburg-Landau theory, which is very similar to the Abclian Higgs model (7). 
Although the details of the dynamics are presumably quite different, the consider- 
ations that led to the flux trapping mechanism (see Section 5.2) are still valid, and 
therefore the mechanism can be tested in superconductor experiments. 

In principle, the required experimental setup is relatively simple. A supercon- 
ducting film is cooled through the phase transition from the normal phase into the 
superconducting phase, and the vortices formed in the transition are detected by 
measuring their magnetic field. 

In practice, there are many difficulties in this kind of an experiment. A slightly 
simpler setup^^ involves a superconducting loop, which is quenched through the 
phase transition. In this kind of a system, the "vortices" show up as a non-zero mag- 
netic flux through the loop. The Kibblc-Zurek mechanism predicts iYxz ~ iL/iY^"^ 
flux quanta, where L is the circumference of the loop, and the magnetic field trapped 
according to the flux trapping mechanism consists of A/trap ~ (e^TL)^/^ flux quanta, 
independently of the quench rate, provided that ^ < L. 

An experiment that was essentially like this was carried out by Carmi et al}'^'^ 
They used a 1 cm^ loop of YBa2Cu307_5 superconductor built from 214 Josephson 
junctions. A light beam was used to heat the loop above ~ 90 K, and when 
the light was switched off, a thermal link to a heat bath cooled the loop through 
the phase transition. The cooling rates varied from 0.3 to 20 K/s. They measured 
the flux through the loop after the transition using a superconducting quantum 
interference device (SQUID) magnetometer placed at a distance of 1 mm from the 
loop, and found that its was zero on average and the distribution had a standard 
deviation of (Xcxp = 7.4 ± 0.7 flux quanta, independently of the cooling rate. 

In this setup, the Kibble-Zurek mechanism predicts that the phase of the order 
parameter is independent in each of the N = 214 segments, and therefore the 
resulting standard deviation should be roughly ctkz = C\fN w 14. 6C, where C is 
some constant of order one. This result agrees well with the observations, and is 
indeed independent of the quench rate, as observed. The contribution predicted by 
the flux trapping mechanism is^'* CTtrap ~ 4 flux quanta, and therefore the experiment 
cannot really distinguish between the two mechanisms. 

What simplifles the superconductor loop experiment is that the "vortices" are 
formed outside the superconductor itself. This, however, also means that the exper- 
iment does not really probe the dynamics inside the superconductor, and indeed in 
the experiment by Carmi et al}^^ the numbers of vortices predicted by the Kibble- 
Zurek and flux trapping mechanisms are both independent of the quench rate. 

Carmi and Polturak^^^ also carried out a similar experiment with a supercon- 
ducting YBa2Cu307_5 film of size 1 cm^. Again, they heated the film with light, 
switched the light off and measured the resulting magnetic flux through the fllm 
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when the system had reached the superconducting phase. They estimated that 
the sensitivity of their measurement was around 20 flux quanta, and they could 
not find any evidence of vortices being formed. They compared this with their 
own theoretical estimate, which was around 10000, and concluded that there is a 
clear discrepancy. However, their theoretical picture was rather different from the 
Kibble-Zurek or flux trapping scenarios, which predict much fewer vortices, of order 
100 and 1, respectively. ^^^'^^ 

In the film experiment, the Kibble-Zurek prediction depends on the quench 
rate, so one could test the scenario by measuring this dependence. However, the 
prediction also depends on details of the dynamics that have not yet been prop- 
erly understood, and therefore it is not a very robust test. Moreover, the flux 
trapping scenario becomes slightly more complicated because the magnetic field ex- 
tends outside the two-dimensional film.^^'^^ For this reason, simple two-dimensional 
simulations such as those in Ref. 64 do not describe the transition correctly. 

On the other hand, as pointed out in Section 5.3, the two scenarios predict very 
characteristic correlations between the produced vortices, and therefore a measure- 
ment of the spatial vortex distribution is the best way of distinguishing between the 
mechanisms. Unfortunately, this was not possible in the experiments carried out 
by Carmi et aZ.^^^'^^° because they only measured the total magnetic flux. 

In order to measure the spatial vortex distribution, one has to face several tech- 
nical problems. First of all, one has to be able to identify individual vortices, which 
means a sensitivity to fluxes at a level well below one flux quantum. Secondly, the 
spatial resolution must be high enough so that one can see the correlations, and 
finally, the measurement must be fast enough so that the vortices do not have time 
to annihilate. There are several different methods available, all of which have 
their strengths and weaknesses. Scanning SQUID microscopy has a very high sen- 
sitivity and reasonably high spatial resolution, but the measurements take minutes. 
Real-time imaging at video rate with a sensitivity to a single flux quantum has been 
achieved using Lorentz microscopy^'"' and magneto-optical imaging. 

The dynamics of the vortices after the transition are complicated by the impu- 
rities in the superconductor fllm, which tend to pin vortices to them. This effect 
can also be used to simplify the measurement by constructing a regular array of 
pinning sites from submicron holes (antidots).^^^ When the vortices get pinned 
to the antidots, their spatial distribution is stabilized and can be measured with 
slower techniques such as scanning SQUID microscopy and scanning Hall probe 
microscopy. 

A further type of an experiment has been proposed by Kavoussanaki et al}'^^ 
who suggested using an annular Josephson junction consisting of two superconduct- 
ing loops separated by a thin layer of insulator. This system is described by the 
sine-Gordon theory, whose fleld is the phase angle difference of the order param- 
eter fields in the two loops. A soliton of the sine-Gordon model corresponds to a 
case with different magnetic fluxes through the two loops, and the Kibble-Zurek 
mechanism predicts that these solitons are formed in the transition. 
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8. Conclusions 

In this article, we have discussed the formation of topological defects in field theory 
phase transitions. In particular, wc have concentrated on the differences between 
gauge field theories and theories with global symmetries. 

While the Kibble-Zurek scenario, ^'^ which applies to theories with global sym- 
metries, has been around for decades and has been tested both in numerical simula- 
tions and in various condensed matter experiments, the same phenomenon in gauge 
field theories has remained poorly understood until very recently. However, in the 
last couple of years, the role of magnetic field has been clarified and a theoretical 
picture of defect formation based on the trapping of magnetic flux has emerged. -'^^ 
Both the Kibble- Zurek and the flux trapping mechanisms are based on a freeze-out 
of long-wavelcngtli modes that are too slow to adapt to the change of the external 
parameters. It is noteworthy that very similar phenomena are expected to take 
place in heavy-ion collisions. 

We have reviewed the basic predictions of the Kibble-Zurek and fiux trapping 
scenarios. The numbers of vortices formed by these two mechanisms depend on 
the temperature, the rate of the phase transition and other parameters in different 
ways. Furthermore, the mechanisms lead to very different spatial distributions of 
vortices: While the Kibble-Zurek mechanism predicts negative correlations between 
the defects, the correlations are positive in the flux trapping scenario. This means 
that flux trapping typically forms clusters of vortices with equal sign, and depending 
on the parameters, these clusters can be large. 

These predictions have been confirmed in numerical simulations, but they can 
also be tested in superconductor experiments, which are made possible by the recent 
developments in experimental techniques. This has been demonstrated in the recent 
pioneering experiments, and there is little doubt that proper experimental 
tests will be carried out in the near future. 

These experiments will provide us with valuable information about the non- 
equilibrium dynamics of gauge field theories during phase transitions. This is ex- 
tremely important for particle physics and cosmology, as well, because once the 
behaviour in a simpler system such as a superconductor is understood, similar the- 
oretical techniques can be applied to heavy ion collisions and cosmological phase 
transitions. 
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